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Dimension 

(L  « length  and  T  a  time,  UP*1) 

AP  Aft  perpendicular 

DBW  Diverging  bow  wave 

DSW  Diverging  stem  wave 

E  Amplitude  of  the  free  wave  spectra 

F  Sine  component  of  the  free  wave  spectra 

Fn  Froude  number,  U/(g  LWL)1/2,  K<,  x  LWL  a  i/Fn2 

FP  Forward  perpendicular 

G  Cosine  component  of  the  free  wave  spectra 


K„  Fundamental  wave  number  (g/U2)  L'1 

LWL  Length  at  waterline  L 

TW  Transverse  wave 

U  Free  stream  velocity  LT"1 

X  Longitudinal  coordinate  (origin  at  FP)  L 

Y  Transverse  coordinate  L 

Z  Vertical  coordinate  L 

g  Gravitational  constant  LT~2 

u  Circular  wave  number  induced  in  y  direction  (sec  9  tan  0) 

w  Circular  wave  number  induced  in  x  direction  (sec  0) 

h  Wave  elevation  L 

q  Direction  of  wave  propagation 

A0  Fundamental  wave  length  (2n/Ko)  L 


0  Angle  from  which  the  transverse  wave  number  is  computed: 

defined  as  the  angle  between  the  direction  of  wave 
propagation  and  the  direction  in  which  the  ship/model 
is  travelling  (see  figure) 
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ABSTRACT 

Wave  patterns  preaicted  by  several  numerical  codes  were  evaluated  by  comparing 
them  with  data  from  model  basin  experiments  on  two  ship  hulls  at  three  Froude  numbers 
each.  In  general,  the  codes  all  underpredict  the  amplitude  of  divergent  waves  springing 
from  the  ships’  bow.  High  wave  number  detail  is  also  lacking  in  the  vicinity  of  the  bow 
wave  cusp  line  compared  to  model  free-sutface  wave  patterns.  Conversely ,  the  codes  tend 
to  overpredict  the  amplitudes  of  waves  close  to  the  ships'  stern  and  in  the  transverse 
wave  system  behind. 

The  codes  are  ranked  according  to  how  closely  they  simulate  the  empirical  results. 
Within  the  ranking,  adjacent  codes  give  similar  wave  predictions  and  the  rank  might  be 
interchanged  for  a  few  Froude  number  cases.  Higher  ranked  codes  gave  consistently  bet - 
ter  predictions  than  the  lower  ranked  codes. 

Problems  for  the  lower  ranked  codes  included:  excessive  wave  damping  such  that 
wa\>es  are  attenuated  near  the  edges  of  the  computational  domain;  reflections  appear  at 
the  outer  boundaries;  high  frequency  “noiss"  exists  beyond  the  nineteen-degree 
envelope  of  the  spreading  wave  train:  and  wave  energy  has  been  severely  underpredicted. 

Standardized  graphic  representations  of  the  numerical  and  experimental  data  are 
included  in  the  appendixes. 

The  predictions  were  made  “blind,"  without  prior  knowledge  of  the  specific 
experimental  results.  In  view  of  the  inherent  difficulties  including  nonlinear  wave  break¬ 
ing  and  real  fluid  effects  ( present  in  the  experiments),  the  agreement  achieved  by  the 
better  codes  is  encouraging. 

ADMINISTRATIVE  INFORMATION 

This  investigation  was  sponsored  by  the  Chief  of  Naval  Research,  Office  of  Naval 
Technology,  Code  ONT^i,  under  the  Surface  Ship  Technology  Program  (ND1A), 

Program  Element  62121N,  Acoustic  and  Wake  Signature  Project  RH21C12,  Task  1, 
Signature  Reduction.  The  DTRC  Work  Unit  numbers  are  1506-710  and  1506-810. 

INTRODUCTION 

The  rapid  increase  in  the  availability  and  power  of  digital  computers  has  made  it 
practical  to  solve  the  Neumann-Kelvin  problem.  Significant  interest  in  the  problem  has 
led  to  three  “International  Workshops"  on  wave  resistance  in  the  last  ten  years.  While  the 
wave  resistance  problem  is  interesting  and  challenging,  the  closely  related  Kelvin  wake 
problem  is  more  exacting  and  has  a  critical  application  in  modeling  ship  wake  electro¬ 
magnetic  signatures,  an  important  practical  problem  since  the  advent  of  satellite  remote 
sensing  of  ship  wakes  which  resulted  from  NASA's  oceanographic  SEASAT  program  in 
1978.  Wave  phase  relations  are  important  regarding  peak  amplitude  and  slope,  hence  the 
occurrence  of  wave  breaking.  Furthermore,  wave  partem  prediction  requires  finer  spatial 
resolution  because  short  wave  lengths  have  an  integral  role  in  radar  cross  section.  Wave 
resistance  prediction  relies  on  an  integrated  quantity  which  is  fairly  insensitive  to  the  high 
wave  number  end  of  the  free  wave  spectrum  nr  phase  relations  in  wave  space. 

Several  researchers  have  developed  numerical  codes  that  may  have  utility  for 
predicting  or  evaluating  the  impact  of  bull  design  on  surface  ship-wake  patterns  regard¬ 
ing  remote  sensing  from  wave-wave  interactions,  wave  breaking,  white-water,  etc.  In 
general,  these  various  codes  have  not  been  validated  with  regard  to  their  overall  power, 
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accuracy,  suitability,  etc.  This  study  was  performed  to  provide  an  burial  validation,  based 
simply  on  experimental  model  basin  surface  ship  wave  pattern  data,  The  objective  is  to 
provide  a  rational  basis  for  selecting  the  code(s)  showing  the  greatest  potential  for  further 
development  as  a  useful  design/evaluation  tool  and  to  set  the  course  and  framework  for 
such  a  development  effort. 

Various  researchers  were  invited  to  submit  numerical  predictions  for  comparative 
evaluation  of  their  codes,  The  specifications  for  the  competition  are  reproduced  in 
Appendix  A.  The  hull  forms  and  Froude  numbers  were  selected  for  specific  cases  that 
had  extensive  Kelvin  wake  data  from  previous  DTRC  model-scale,  towing  tank  experi¬ 
ments,  These  data  were  not  disclosed  to  the  researchers;  thus  they  were  afforded  a  “blind” 
test  for  their  code's  predictive  capability.  Three  Froude  numbers  for  each  of  two  nonpro- 
pelled  ship  models,  DTRC  Model  3531  (QUAPAW)  and  Model  5415,  were  selected  for 
the  “competition.”  Most  invited  entries  considered  only  the  QUAPAW  model  because 
they  were  unable  to  deal  with  the  transom  stem  of  Model  5415,  a  high  speed  destroyer- 
type  hull.  Standard  formats  were  established  for  presenting  the  codes’  predictions  (see 
Appendix  A)  so  that  comparisons  might  be  facilitated  by  overlaying  transparencies  based 
on  the  experimental  results  presented  m  Appendix  B. 

Three  “referees”  evaluated  the  entries  independently  and  then  met  to  achieve  a  con¬ 
sensus  regarding  the  weaknesses  and  strengths  of  each  entry.  Areas  of  consideration 
included:  wake  realism,  relative  wave  amplitudes,  appropriate  wave  radiation  and/or 
attenuation,  and  free-wave  spectra. 

The  following  sections  discuss  the  ship  hull  characteristics,  the  analytical  models, 
rationales  for  the  data  formats  and  comparative  evaluation,  results,  and  conclusions.  The 
graphic  empirical  data  and  each  numerical  code  entry  are  presented  in  individual  appen¬ 
dixes,  in  the  order  of  their  assessed  proficiency. 

PHYSICAL  MODEL  CHARACTERISTICS 

The  QUAPAW  (ATF  110)  is  an  ocean-going  fleet  mg  represented  by  the  l/12th 
scale  model,  DTRC  Model  3531,  appended  with  a  rudder.  The  model’s  wa'erline  length  is 
16.25  ft,  so  the  three  Froude  numbers,  0.2131,  0.25,  and  0.3197  correspond  to  full-scale 
ship  speeds  of  10.0, 11.7,  and  15.0  knots,  respectively.  The  top  ship  speed  is  about  16 
knots.  This  ship  is  of  particular  interest  because  full-scale  data  from  the  so-called 
“Georgia  Straits”  experiment,1  are  available  for  further  validation. 

DTRC  Model  5415  represents  a  prototype  twin  screw,  transom  stem  destroyer.  This 
model  is  fitted  with  a  large  sonar  dome.  The  waterlme  length  of  this  model  is  18.767  ft; 
thus  the  Froude  numbers,  0.25,  0.2755,  and  0.4136  correspond  to  ship  speeds  of  18.1, 
20.0,  and  30.0  knots,  based  on  a  24.822  scale  ratio.  This  model  was  appended  with 
shafts,  struts,  and  rudders. 

The  model  experiments  were  conducted  with  freedom  in  pitch  and  heave  so  that  the 
hull  took  on  a  trim  condition  when  underway.  The  trim  data  were  supplied  to  the 
numerical  modellers  so  that  wave  patterns  would  be  derived  for  the  "as  trimmed  condi¬ 
tion."  Similarly,  the  “lines”  for  the  two  models  were  furnished  in  the  form  of  ship 
panelization  files,  depicted  graphically  in  Figs.  1  and  2  for  the  two  ship  model  hulls.  The 
numerical  modellers  were  not  restricted  to  a  specific  panelization  scheme,  however,  and 
some  used  modified  panels. 
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Fig.  1.  Typ<cal  pane!  scheme  for  Model  353i . 


Fig.  2,  Typical  panel  scheme  lor  Model  5415, 


ANALYTICAL  MODELS 


The  predictive  programs  which  were  evaluated  can  be  divided  into  two  categories: 
those  using  Havelock  sources,  which  automatically  satisfy  the  linearized  free-surface 
condition,  distributed  over  the  surface  of  the  body;  and  those  which  use  Rankine  sources 
distributed  over  the  body  and  the  undisturbed  ffee-surface.  These  latter  methods  are 
closely  related  to  Dawson's  pioneering  approach2  wherein  singularities  are  distributed  on 
the  free  surface  (as  well  as  on  the  hull  surface)  in  such  a  fashion  that  they  satisfy  a 
linearized  free  surface  condition. 

The  creators  of  five  computer  programs  based  on  the  Havelock  source  potential 
were  asked  to  make  Kelvin  wake  predictions.  These  authors,  their  employers,  and  their 
computer  programs'  names  are  given  in  Table  1.  Creators  of  seven  Rankine  Source  pro¬ 
grams  were  aiso  asked  to  make  Kelvin  wake  predictions.  A  list  of  these  authors,  their 
employers,  and  the  names  of  their  computer  programs  are  given  in  Table  2. 

Table  1.  Computer  programs  based  on  Havelock  source  potential. 


I  Author 

Employer 

Program  Name 

R.F.  Beck 

University  of  Michigan 

DOCTORS  N-  K3'4 

F.  Noblesse 

DTRC 

EPPAC1,  FARWAV5 

T.  Schmidt,  J.V.  Rattayya, 

Lockheed,  Advanced 

and  J.R.  Brooks 

Marine  Systems 

C.A.  Scragg 

SAiC,  San  Diego 

NK-  2* 

J.  Telste 

DTRC 

SWIM,  SWIMFS7 

Table  2.  Computer  programs  based  on  Rankine  source  potential. 


Author 

Employer 

Program  Name  i 

W.  Cheng 

DTRC 

XYZFS8 

J.  Dean  (proprietary) 

DTRC 

— 

Y.H.  Kim  and  S.H.  Kim 

DTRC 

SWIFT* 

L.  larsson 

SSPA 

—  - 

B.  Rosen 

Self 

FLOPAN10 

P.  Sclavounos  and 

J.N,  Newman 

MIT 

NKSHIP11 

K.  Weems  and  C.  Oliver 

SAIC, 

Annapolis 

—  - 

On  12  and  13  January  1988,  a  workshop  was  held  at  DTRC  to  discuss  and  comoare 
results  of  the  Kelvin  wake  computations.  A  formal  letter  announcing  the  workshop  was 
sent  to  each  participant  The  letter  is  reproduced  in  Appendix  A.  This  workshop  was  at¬ 
tended  by  the  various  individuals  who  had  made  Kelvin  wake  computations,  and 
additional  corporate  representatives  from  many  of  the  organizations  represented  in  the 
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computational  effort,  la  addition,  the  workshop  was  held  in  conjunction  with  a  meeting 
of  the  Society  of  Naval  Architects  and  Marine  Engineers  (SHAME)  Analytical  ShipWave 
Relations  (H-5)  panel.  The  participation  of  the  H-5  panel  neatly  broadened  the  atten¬ 
dance  at  the  meeting,  both  in  numbers  of  people  and  in  technical  breadth, 

During  the  course  of  the  workshop,  the  authors  of  the  various  codes  discussed  then- 
prediction  methods  and  our  comparative  evaluations  of  the  codes  were  presented. 
Graphic  representations  of  the  various  invited  predictions  are  reproduced  in  Appendixes 
C  -  K  along  with  written  discussions  that  the  authors  prepared  for  the  workshop  in  re¬ 
sponse  to  the  “Wake-Off  Questionnaire”  and  “Kelvin  Wave  Computations"  (reproduced 
in  Appendix  A.) 

Although  it  was  hoped  that  Dr.  Lars  Larsson  of  SSPA  and  Chalmers  University 
would  be  able  to  contribute  to  the  computational  effort,  his  travel  commitments  made  it 
impossible  for  him  to  complete  the  computations  within  the  time  frame  required.  In  addi¬ 
tion,  two  of  the  programs  (those  of  Schmidt,  et  al,  and  Weems  and  Oliver)  were  judged  to 
need  further  development,  and  thus  these  computer  programs  are  not  included  in  the 
comparative  evaluation. 


DATA  FORMATS 

Details  of  the  formats  that  were  specified  for  graphical  data  presentations  are  shown 
in  Appendix  A.  Three  types  of  plots  were  included;  iso-elevation  contours,  longitudinal 
wave  cuts,  and  free-wave  spectra.  Scales  were  selected  based  on  the  model  test  data 
ranges  and  more-or-less  standard  dimensionless  parameters. 

The  most  complete  spatial  representation  of  free-surface  elevation  are  contour  plots 
showing  lines  of  iso-elevation  at  a  specific  interval;  c.g.,  0.01/Ao,  where  K0  is  the 
fundamental  wave  number  and  l/K0  =  2h0  or  twice  the  pressure  head  of  the  flow).  The 
experimental  contours,  obtained  from  stereophotographic  data,  were  all  drawn  with  a  2% 
(AZK0  -  0.02)  interval  as  a  compromise  between  too  many  confusing  lines  (smaller 
interval)  and  insufficient  detail  (larger  interval), 

Basic  two-dimensional  wave  profiles,  obtained  with  capacitance  wave  probes,  are 
depicted  m  longitudinal  cuts  of  the  free  surface  elevation  in  a  vertical  plane  parallel  to  the 
ship’s  longitudinal  centerline.  The  transverse  location  for  wave  cuts  is  Y  *  0.324*LWL 
for  Model  5415  and  Y  =  0.373*LWL  for  Mode!  3531.  The  disturbed  wave  elevation  (Z) 
and  longitudinal  coordinate  (X)  have  been  normalized  by  the  fundamental  wave  number, 
i.e.,  Z*K0  and  X*K0.  The  origin  is  abreast  the  ship's  forward  perpendicular,  at  the  undis¬ 
turbed  free-surface,  as  sketched  below.  Note,  the  vertical  scale  is  greatly  magnified 
compared  to  the  horizontal  scale,  so  that  the  wave  profiles  have  a  very  steep  appearance. 
The  range  for  the  vertical  scale  was  specified  to  be  ±  0.1  %,  or  ±  20  %  of  the  fluid 
velocity  head;  experimental  wave  cut  surface  elevations  were  bounded  by  this  range. 
The  horizontal  scale  ranged  from  0  to  100,  so  that  roughly  sixteen  fundamental  wave¬ 
lengths,  are  contained  in  the  graph  Q.0  =  2ji/K0). 

In  order  to  make  a  quantitative  comparison  between  analytical  predictions  and  ex- 
peiimental  data  it  is  important  to  quantify  the  uncertainiues  generated  during  the  model 
experiment.  Although  there  are  a  vast  number  of  potenual  sources  of  inaccuracy  only  a 
few  of  the  more  important  will  be  quantified  here;  model  speed,  model  construction, 
model  surface  finish,  and  the  repeatability  of  the  capacitance  probe  data. 
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The  speed  of  the  model  being  towed  on  the  carriage  is  measured  using  a  toothed 
gear,  which  rotates  along  with  one  of  the  carriage  drive  wheels  and  triggers  a  magnetic 
pick-up  with  a  pulsed  output.  The  resultant  precision  error  of  the  measurement  of  car¬ 
riage  speed  is  ±  0.001  ft/sec  and  the  bias  error  associated  with  that  is  0.005  ft/sec. 

Botn  the  models  which  were  used  as  part  of  the  this  study  were  of  wooden  construc¬ 
tion  which  ensures  a  certain  stability  of  shape  over  time.  The  model  section  shape 
contours  are  constant  to  within  ±  1/16  in.  and  the  “waviness"  is  less  that  1/32  in.  between 
ends  of  a  3-ft  batten.  Model  5415  was  a  newly  constructed  model,  and  thus  the  above 
tolerances  apply,  however  Model  3531  was  built  in  the  1940’s  and  most  probably  has 
additional  distortion  enough  to  double  the  tolerances  in  section  shape  contours  from  that 
of  a  new  model. The  surface  finish  for  both  models  is  obtained  using  enamel  paint  applied 
to  a  high  gloss  smooth  finish  with  roughness  over  the  surface  of  no  more  than  60-100 
micro-inches  RMS. 

In  an  effort  to  examine  the  repeatability  of  the  capacitance  probe  response  to 
waves  produced  by  the  model  during  a  run,  two  repeat  wave  records  from  the  same  probe 
were  compared.  For  two  model  speeds,  4.36  knots  and  11.65  knots,  these  comparisons 
yielded  a  relative  precision  error  (standard  deviation)  of  2  percent  and  1.8  percent,  re¬ 
spectively,  referenced  to  the  peak  amplitude  of  the  wave  record.This  error  was  obtained 
by  calculating  a  standard  deviation  between  two  wave  records  using  each  of  the  2000 
data  points  in  the  record.  A  more  detailed  account  of  the  uncertainties  introduced  in  ob¬ 
taining  measurements  of  model  wave  heights  is  found  in  Reference  12. 

Free  wave  spectra  show  the  distribution  of  wave  energy  density  amplitude  and  its 
two  (Fourier)  components  as  a  function  of  the  transverse  wave  number,  u.  Transverse 
wave  energy  is  manifested  in  the  very  low  wave  number  range  and  divergent  waves  cor¬ 
respond  to  the  remaining  wave  numbers.  Note,  the  three  graphic  representations  for  each 
data  set  are.  in  some  sense,  redundant— one  can,  in  theory,  derive  each  graph  from  either 
one  of  the  other  two. 


6 


EVALUATION  PROCESS 

The  enfriei  were  evaluated  independently  by  three  referees  (the  authors), 
against  standardized  wake  criteria  agreed  upon  at  the  beginning  of  the  evaluation  pro* 
cess.  They  then  met  to  compare  then  individual  findings  and  achieve  a  consensus 
regarding  the  weakness  and  strength  of  each  entry.  Areas  of  consideration  included: 
wake  pattern  realism  (appropriate  cusp  wave  angle,  phase  relationships,  fine  scale  detail, 
and  no  unrealistic  reflection  planes):  relative  amplitudes  of  bow  divergent,  stem  divergent 
and  transverse  waves;  appropriate  wave  radiation  trajectory  and  attenuation,  and  free 
wave  specna  distributions  with  regard  to  amplitude  and  phase. 

The  findings  are  discussed  in  the  .  uwing  sections  with  tables  of  comparative 
wave  amplitude  and  relative  ranking.  For  completeness,  the  entries  are  all  reproduced  in 
appendixes,  in  the  order  of  their  closeness  of  fit  to  the  experimental  wave  pattern  wakes. 
For  this  assessment,  no  account  was  taken  of  the  relative  software  documentation  and 
user  friendliness,  nor  of  the  hardware  and  CPU  time. 

RESULTS 

The  evaluations  were  based  largely  on  graphic  plot  comparisons.  The  reader  is 
invited  to  make  his  own  set  of  clear  “viewgraphs*  from  the  experimental  data  represen¬ 
tations  in  Appendix  B.  These  can  then  be  overlaid  on  corresponding  numerical  data  plots 
to  corroborate  the  results  discussed  below. 

GENERAL  OBSERVATIONS 

There  were  some  basic  similarities  between  the  results  Of  many  of  the  computer  pro¬ 
grams.  In  general,  the  agreement  between  prediction  and  measurement  was  better  at  the 
highest  Froude  number  for  each  model.  Other  similarities  were  found  in  the  bow  and 
stem  wave  systems:  the  amplitude  of  the  bow  wave  was  consistently  underpredicted  and 
the  amplitudes  of  the  stem  and  transverse  waves  were  overpredicted.  Table  3  shows 
several  comparisons,  derived  from  wave  cut  graphs,  for  predicted  wave  amplitude 
divided  by  the  corresponding  experimental  wave  amplitude.  This  ratio  is  seen  to  be  less 
than  one  for  the  peak  value  arising  from  the  diverging  bow  wave  (DBW)  and  greater 
than  one  for  most  peak  diverging  stem  wave  (DSW)  predictions.  The  amplitude 
mismatch  for  the  predictions  made  on  diverging  stem  waves  (and  transverse  wave  trains) 
is  decidedly  worse  for  Model  3531  than  for  Model  5415. 

The  crest  lines  along  the  bow  wave  cusp  line  are  predicted  to  have  larger  angles 
than  shown  experimentally.  This  may  be  caused  by  the  transverse  waves  being  overpre- 
dicted,  while  DBW  is  underpredicted. 
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Tibie  5.  Relative  wave  amplitudes. 


TabieSa.  Model  5415. 


CODE 

NAME 

0.25 

FROUDE  NUMBER 

0.2755 

0.41360 

DBW 

DSW 

TW 

DBW 

DSW 

TW 

DBW 

DSW  TW 

FARWAV 

0.7 

1.36 

2,0 

0.6 

1.45 

3.5 

0.8 

1.0  1.4 

NK-2 

0.2 

1.45 

3.0 

0.33 

1.80 

3.0 

0.67 

1.26  1.4 

SWIFT 

0.2 

0,64 

— 

0.25 

0.62 

0.45 

0.80  — 

XYZFS 

0.16 

0.05 

— 

0.13 

0,06 

~~ 

0.45 

0.73 

Note:  Values  show  the  ratio  of  predicted  wave  amplitude  to  experimental  wave  ampli¬ 
tude.  DBW:  diverging  bow  wave,  DSW:  diverging  stern  wave,  TW:  transverse  wave. 


Table  3b.  QUAPAW,  Model  3531. 


CODE 

FROUDE  NUMBER 

NAME 

0.2131 

0.25 

0.3197 

DBW 

DSW 

TW 

DBW 

DSW 

TW 

DBW 

DSW 

TW 

FARWAV 

0.86 

2.7 

4.0 

0.80 

6.2 

5.0 

1.0 

1.8 

2.0 

NK-2 

0.52 

2.6 

2.5 

0.6 

10 

7.5 

0.75 

2.4 

3.5 

SWIFT 

0.25 

0.92 

— 

0.37 

1.8 

— 

0.56 

1.33 

— 

FLOPAN 

0.36 

0.35 

— 

0.47 

11 

— 

0.71 

0.89 

■ — 

XYZFS 

0.14 

0.16 

0.05 

0.45 

0.88 

2.5 

0.32 

1.26 

1.8 

Dean 

0.11 

0.05 

0.20 

0.18 

0.16 

1.7 

0.38 

0.44 

2.5 

SWIMFS 

0.93 

1.8 

2.4 

0.58 

6.2 

8 

— 

4 

7 

Note:  Values  show  the  ratio  of  predicted  wave  amplitude  to  experimental  wave,  amplitude. 
DBW:  diverging  bow  wave,  DSW:  diverging  stern  wave,  TW:  transverse  wave. 

A  strong  intermediate  crest  line  is  seen  in  the  empirical  wave  patterns  in  the  region 
between  the  diverging  bow  and  stem  wave  systems.  It  may  be  emanating  from  the  for¬ 
ward  shoulder  region  of  the  hull.  This  feature  is  generally  not  seen  in  predicted  wave 
patterns;  it  is  an  important  feature  because  the  local  wave  slope  (and  thus  wave  breaking) 
is  affected  by  such  pattern  details.  This  feature  is  evident  in  the  experimental  wave  cuts  as 
an  “extra"  wave  between  the  bow'  and  stem  divergent  wave  peaks. 

/mother  weakness  evident  in  the  code  predictions  (except  for  FARWAV)  was  that 
wave  amplitude  is  over-attenuated  downstream,  along  the  bow  wave  cusp  line.  The 
Rank  me  source  programs  generally  had  greater  damping  of  waves  diverging  away  from 
the  ship.  At  the  extreme,  FLOPAN  predictions  showed  no  waves  at  the  boundaries  of  the 
computational  domain.  Conversely,  experimental  transverse  waves  were  attenuated  more 
rapidly  than  predictions  (and/or  the  classical  analytical  R“1/:  dependence). 
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Experimental  free-wave  amplitude  spectra  exhibit  more  sharply  defined  nodes  or 
null  points;  hence,  wave  phase  relationships  were  incorrectly  predicted  numerically. 

SPECIFIC  COMPARISONS  AND  RANKINGS 

Comparative  rankings  of  the  codes  are  enumerated  in  Table  4  for  three  categories: 
contour  plots,  wave  cuts,  and  free-wave  spectra.  Some  entires  did  not  include  one  catego¬ 
ry  of  graph  (e.g.,  spectra);  they  were  judged  solely  by  the  categories  submitted  (wave  cuts 
and  contour  plots),  without  prejudice. 

Table  4.  Code  ranking  by  category. 


Table  4a.  Model  5415. 


Contour  Plot 

Wave  Cut 

Spectra 

1.  FARWAV 

1.  FARWAV 

1.  FARWAV 

2.  NK-2 

2.  SWIFT 

2,  NK-2 

3.  SWIFT 

3.  NK-2 

3.  XYZFS 

4,  XYZFS 

4.  XYZFS 

1 

Table  4b.  QUAPAW,  Model  3531. 


Contour  Plot  .  Wave  Cut  Spectra 

1.  FAR  WAV,  FLOPAN*  1.  FAR  WAV  iT  FAR. WAV 

2.  SWIFT,  NKSHEP,  2.  NK-2  2.  DOCTORS  N-K 

NK-2, DOCTORS  N-K 

3.  Dean,  SWIM,  XYZFS  3.  DOCTORS  N-K  3.  NK-2 

4.  SWIFT,  NKSHIP  4 .  SWIM,  XYZFS 

5.  XYZFS  5.  Dean 

■  6.  FLOPAN,  Dean. 

!  SWIM 

|  *Note:  FLOPAN  predictions  are  supenor  very  close  to  the  hull,  while  FARWAV 

| _ predictions  are  better  for  waves  radiating  away  from  the  hull.  _ 

Five  of  the  computational  results  were  significantly  better  than  the  others,  those  of: 
FARWAV.  NK-2.  DOCTORS  N-K.  SWIFT,  and  FLOPAN.  FARWAV  predictions  were 
consistently  superior  for  both  models  at  all  speeds  and  iB  the  three  evaluation  categories. 
Wave  cuts  were  very  similar  between  FARWAV  and  NK-2  owing  to  similar  mathematical 
representations  of  source  potentials.  However,  spectra  predicted  by  FARWAV  compar  ed 
more  favorably  to  the  experimentally  determined  resuits;  the  diverging  waves  did  not 
decay  as  fast  as  those  from  the  other  Havelock  source  methods. 

It  is  interesting  to  note  that  FARWAV  used  a  more  rudimentary  slender-ship 
(zeroeth-order)  approximation  for  the  source  distributions  on  the  ship  hull.  Subsequent 
to  the  Wake-Off  workshop,  NK-2  was  exercised  with  the  slender-ship  approximation 
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and  the  results,  which  essentially  duplicated  the  FARWAV  predictions,  are  included  in 
Appendix  L. 

As  stated  previously,  PLOPAN  wave  patterns  died  out  much  too  rapidly,  However, 
FLOPAN  gave  the  best  predictions  of  all  for  the  details  of  the  free-surfaee  disturbance  in 
the  region  of  one  ship’s  beam  from  the  model.  These  predictions  showed  frcc-surface 
details  (evident  in  experiments)  which  were,  at  best,  hinted  at  by  the  predictions  from  the 
better  of  the  other  programs.  Subsequent  to  the  workshop,  FLOPAN  was  rerun  with  two 
alternative  finite  differencing  schemes  that  served  to  alleviate  the  wave  attenuation  prob¬ 
lem  and  enhanced  the  detail,  but  with  some  signs  of  instability.  These  modified  results  are 
included  in  Appendix  M. 

The  predictions  by  SWIFT  were  very  similar  to  FLOPAN,  with  high  attenuation  at 
the  boundaries  of  the  computational  domain.  The  wave  patterns  near  the  ship  lacked 
some  of  the  fine  scale  detail  seen  in  FLOPAN’s  contour  plots.  NKSHIP  predictions  seem 
qualitatively  similar  to  SWIFT’s,  but  only  a  limited  range  was  submitted  so  our  finding  is 
somewhat  tentative  in  this  case. 

The  remaining  codes  gave  noticeably  poorer  quality  results.  XYZFS  predictions  for 
both  models  are  marginally  acceptable  at  the  highest  Froude  numbers.  XYZFS  was  rerun 
after  the  workshop  with  a  revised  surface  panel  distribution.  Improved  results  are  very- 
similar  to  those  from  SWIFT,  and  are  included  in  Appendix  N. 

Dean’s  results  with  the  proprietary  code  exhibited  unnatural  growth  in  transverse 
waves  downstream  at  low  Froude  numbers.  SWIM/SWIMFS  predictions  appear  to  have 
problems  with  numerical  resolution  and  waves  propagating  upstream. 

DISCUSSION 

Havelock  source  codes  generally  Out  performed  Dawson/Ranlrine  source  codes.  On 
the  other  hand,  near-field  results  such  as  FLOPAN’s  may  favor  the  Rankine  source 
approach. 

Wave  profiles  along  the  hull  surface  are  readily  available  and  would  have  been 
useful  for  comparing  “very”  near-field  predictions.  This  output  was  overlooked  when 
setting  the  ground  rules  for  the  competition  (empirical  wave  profiles  were  unavailable  as 
ground  truth,  however).  Similarly,  predicted  wave  pattern  resistance,  Cw.  should  have 
been  requested  in  the  specification.  These  data  are  shown  in  Table  5  for  those  cases  where 
submitted  by  the  invitees. 
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Tibi#  5,  Wav#  pattern  resistance  coefficient  predictions.  (Values  shown  are  0w  x  1000) 


Model  5415 

Model  3531 

Code 

Fn  - 

0.25 

0.2755  0.4136 

0.2131 

0.25 

0.3197 

DOCTORS  N-K 

0.96 

3.94 

6,0 

FLOPAN 

0.21 

0.62 

3.5 

NK-2 

0.51 

1.39  4.4 

1.37 

6.40 

7.9 

XYZFS 

0.80 

0.89  2.3 

0.66 

0.48 

5.1 

Experiment 

0.37 

0.46  2.4 

0.21 

0.40 

1.8 

Some  of  differences  between  experimental  and  predicted  wave  patterns  can  be  at¬ 
tributed  to  real  fluid  or  viscous  effects.  The  experimental  stem  wave  cusp  line  is 
noticeably  offset  outboard  of  the  model  compared  to  the  predicted  wave  patents.  Some  of 
differences  between  experimental  and  predicted  wave  patterns  can  be  attributed  to  real 
fluid  or  viscous  effects.The  experunetnal  stem  wave  cusp  line  is  noticeably  offset  out¬ 
board  of  the  model  compared  to  the  predicted  wave  patterns.  The  thick  boundary  layer  at 
the  stem  is  probably  affecting  the  pattern  here  as  well  as  contributing  to  the  reduction  in 
stem  wave  amplitudes. 

Attenuation  of  the  transverse  waves  may  likewise  be  affected  by  the  viscous  wake 
trailing  behind  the  experimental  model.  In  this  regard,  viscous  effects  are  exaggerated  (by 
the  low  Reynolds  number)  in  the  experiments  compared  to  full-scale  ships. 

Science  Applications  International  Corp.  (SAIC)  has  performed  detailed  analysis  of 
full-scale  wake  data  for  the  ATF  QUAPAW.13  Figure  3  shows  SAlC's  (ship)  free-wave 
amplitude  spectrum  compare!  with  Model  3531  data  for  an  unpropelled  15-knot  case  (Fn 
*  0.32).  The  full-scale  data  represent  a  15.5-knot  run;  specific  displacement  and  trim  are 
unknown,  but  believed  to  be  relatively  near  those  used  in  the  experiment.  In  view  of  the 
difficulties  inherent  in  full-scale  experimentation,  the  agreement  seen  for  these  two  data 
sets  is  encouraging. 

The  codes  generally  assume  that  the  ship  hull  is  “wall  sided"  at  the  free-surface, 
i.e.,  the  verocal  gradient  in  the  hull  offsets  is  zero.  This  simplification  may  affect  the 
predicted  wave  patterns. 

CONCLUSIONS 

One  might  have  hoped  for  better  results  from  the  comparison  of  the  experimental 
and  computational  results.  The  technical  problems  are  difficult  and  wide  ranging,  but  are 
not  intractable.  The  work  of  Noblesse,  Scragg,  Beck,  Kim,  and  Rosen  indicate  that  there 
is  a  good  chance  for  successfully  predicting  Kelvin  wave  fields,  particularly  if  a  method 
can  be  developed  which  can  take  advantage  of  the  excellent  near-field  calculations  of 
FLOPAN  and  the  superior  faMield  capability  of  the  Havelock  source  methods.  The 
Wake-Off  workshop  helped  to  focus  on  the  directions  for  future  improvements. 

The  fact  that  Noblesse,  Scragg,  and  Beck  can  produce  results  which  are  virtually 
identical  when  they  use  consistent  approximations  of  the  body  boundary  condition  is  cer¬ 
tainly  an  encouraging  finding.  This  is  particularly  true  when  one  considers  the  large 
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amount  of  scatter  in  wave  resistancecalcnlations  which  were  performed  as  recently  as 
four  or  five  years  ago,  It  is  to  be  expected  that  Noblesse,  Scragg,  and  Beck  would  all  pre¬ 
dict  similar  wave  resistance  values;  this  implies  an  unprecedented  degree  of  consistency 
in  their  solution  to  the  Neumann-Relvin  problem,  Although  it  should  not  be  given  too 
much  significance,  it  is  also  encouraging  that  the  phases  of  the  humps  and  hollows  of  the 
free=wave  spectra  predicted  by  Noblesse,  Scragg,  and  Beck  are  essentially  in  agreement 
with  the  experimental  data,  although  the  amplitude  does  not  always  agree.  Less  satisfac¬ 
tory  results  are  obtained  when  the  sine  and  cosine  components  of  the  free-wave  spectra 
are  compared  with  the  model  data. 

Seragg’s  studies  of  the  effects  of  smaller  panels  near  the  waterline  show  a  disturbing 
divergence  of  the  wave  resistance  due  to  the  behavior  uf  the  waterline  integral.  This  is 
caused  by  the  fact  that  the  source  strengths  on  the  panels  near  the  free  surface  are  ap¬ 
proaching  infinity  as  the  panel  size  decreases;  so,  while  the  product  of  the  source  strength 
and  the  panel  area  remains  bounded,  the  source  strength  in  the  waterline  integral  goes  to 
infinity.  This  leads  to  anomalous  results  for  the  wave  resirance  and,  obviously,  for  other 
quantities  that  may  be  less  sensitive  to  the  waterline  integral  than  is  the  wave  resistance. 

A  close  examination  of  the  derivation  of  the  waterline  integral  is  probably  warranted,  as 
is  exploring  ways  to  evaluate  the  waterline  integral  without  using  the  source  strength  on 
the  free  surface.  Perhaps  a  method  based  purely  on  the  wave  potential  would  avoid  these 
difficulties. 

Additional  work  is  still  required  to  obtain  an  efficient,  fast,  and  accurate  evaluation 
of  the  Green’s  function  of  the  Havelock  source  potential.  Other  technical  areas  needing 
further  effort  include:  the  proper  treatment  of  non-wall-sidedness,  the  dependence  of  the 
solution  on  panelization  or  convergence,  the  effect  of  propulsion  on  the  wave  system,  the 
proper  boundary  conditions  to  apply  with  a  transom  stem,  and  the  inclusion  of  lifting 
components.  An  ultimate  task  would  be  the  inclusion  of  the  full  nonlinear  free-surface 
condition,  at  least  near  the  ship. 


Fig.  3.  Measured  wave  spectra  tor  the  ATF  QUAPAW  and  Model  3531 . 
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APPENDIX  A 

SPECIFICATIONS  FOR  WAKE-OFF  COMPETITION 

TASK 

Use  your  free-surface  Kelvin  Wake  code  to  compute  the  free  surface  elevations  in 
the  near  field  of  two  models,  one  a  high  speed  transom  stem  destroyer  (Modei  5415)  and 
the  other  the  QUAFAW  (Model  3531),  The  calculations  should  produce  the  entire  wave 
system  in  the  area  bounded  by  a  cross-track  line  0.1  ship  iength  ahead  of  the  bow  and 
another  ctoss-  track  line  three  ship  lengths  downstream  from  the  stem.  The  calculations 
should  be  made  for  several  values  of  Froude  number  for  each  model,  as  specified  in  Table 
A.l.  The  calculations  should  be  presented  as  surface  contours,  as  wave  cuts  parallel  to 
the  ship  track,  and  as  spectral  representations  of  the  wave  cuts.  The  calculations  should 
be  presented  for  the  parametric  values  shown  in  Fig.  A.l. 

Table  A.l,  Model  experimental  condl'rins. 


Model  5415 

LWL  =  18.787  ft 

X  m  24.8220 

V*  (knots) 

Fh 

Depth  to  Keel  (ft) 
at  FP  at  AP 

0.0 

1.0000 

20.22 

20.22 

18.1 

0.2500 

21.00 

20.30 

20.0 

0.2755 

21.20 

20.70 

30.0 

0.4136 

20.00 

23.60 

Wave  cut  at  y/LWL  2= 

6.083/18.767  2  0.324 

Model  3531  (QUAPAW) 

LWL  *  16.250  ft 

A  ■  12.000 

Vk  (knots) 

Fr. 

Depth  to  Keel  (ft) 
at  FP  at  AP 

0.0 

0.0000 

13.60 

16.75 

10.0 

0.2131 

13.87 

17.28 

11.7 

0.2500 

13.94 

17.52 

15.0 

0.3197 

14.43 

18.42 

Wave  cut  at  y/LWL  = 

6.060/16.250  =  0.373 
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Fig  A.lb.  Free  wave  spectra 


1.  (Continued) 


OBJECTIVE 

The  results  will  be  evaluated  with  a  view  toward  using  the  code  in  the  development 
of  the  Surface  Ship  Hydrodynamic  Wake  Reduction  Code.  The  evaluation  will  include 
the  following  criteria: 

1.  For  Accuracy, 

Wave  length 

Elevation  pattern 

Wave  cut  spec  ha 

Wave  amplitude  (Incal  slope) 

Peak  amplitude 

Minimum  resolvable  wave  length 

2.  For  Physical  Code 

User  friendliness 
Overall  power 

Efficiency  — *  fast  and  easily  modified  input  preparation 
Flexibility  for  modifying  geometry 
Modular  format 

Adaptability  to  different  computers 
Ability  to  provide  a  graphical  display 
Ease  in  obtaining  arbitrary  wave  cut  spectra 
Accommodation  for  nonlinear  effects 
Wide  range  of  Froude  numbers 

3.  For  Documentation 

Availability 

Clarity 


GEOMETRY 

The  geometry  for  paneling  the  hull  is  available;  calculations  should  be  done  for  the 
sinkage  and  irun  conditions  found  in  Table  A.l. 

FORMAT 

The  results  are  to  be  presented  in  graphical  format  as  shown  in  Fig.  A.l  so  that  com¬ 
parisons  between  calculated  and  empirical  data  may  be  made  easily  and  directly  by 
overlaying  the  graphics.  The  data  are  to  be  normalized  in  the  manner  shown  in  Fig.  A.l. 
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DEPARTMENT  OF  THE  NAVY 

BAVie  W.  TAVLOrt  HAVAL  SHIP  PiSCAHeM 
AH  0  Di Vi Lo PUSWT  eUHTtm 
HiAbeUAftfSiti 

•  CtHEibA,  MARYLAND  adeiA^floo 


ANNAP0U*  UAlbMAYSilV 
AnNaPoUI,  Ud  114dM64t 

CAADIAaCK  LAiOH  At  o  H  V 
SEYWSSoA.  Mo  ||944Heoa 


IN  NBpLY  1£?Sh  Yoi 

9  Dec  1987 


Dear 

As  you  are  aware  from  previous  ccrrmunleatlons ,  we  are  planning  to  hold  a 
two-day  workshop  to  discuss  the  results  of  the  comparative  Kelvin  wake  calcula¬ 
tions  on  12  and  13  January  1988.  You  are  invited  to  attend  and  participate  in 
this  workshop,  which  will  also  be  a  joint  meetina  with  the  Sh'AME  H-S  ANALYTICAL 
SHIP  WAVE  RELATIONS  PANEL. 


The  exact  agenda  of  the  meeting  has  not  yet  been 
general  outline  Is  to  have  these  participants  who  hav 
either  the  Cwapaw  or  DDG-S1  make  approximately  a  30  m 
details  cf  their  computational  method.  It  is  to  be  a 
are  familiar  with  the  boundary  value  problem  assoclat 
problem  for  a  ship  in  steaay  motion  on  a  free  surface 
and  Rankine  source  methods  for  solving  this  problem, 
niques  and  results  will  then  be  presented.  This  will 
tations  on  what  was  learned  in  doing  the  computations 
difficulties  which  were  encountered.  This  presentati 
those  who  did  calculations.  Finally,  the  comparative 
putational  results  will  be  presented  followed  by  a  se 


prepared,  however,  the 
e  prepared  results  for 
irute  presentation  on  the 
ssumed  that  all  participants 
ed  with  the  Kelvin  wake 
,  and  the  Havelock  source 
The  experimental  tech- 
be  followed  by  presen- 
and  the  specific 
or,  will  be  optional  for 
evaluations  of  the  com- 
t  of  general  conclusions. 


For  these  individuals  who  have  done  comparative  computations,  It  Is  asked 
that  you  present  two  or  three  vieworaphs  at  the  beginning  of  your  presentation 
which  answer  the  questions  that  are  attached  as  Enclosure  (1).  The  purpose  of 
these  questions  Is  to  provide  a  uniform  set  of  measures  as  to  the  maturity  and 
sophistication  of  each  of  the  codes.  In  your  presentation  of  the  details  of 
your  code,  It  Is  suggested  that  you  Include  Information  On  the  Green's  function 
used,  the  Integration  method  over  the  panels,  the  method  for  satisfying  the  body 
and  free  surface  boundary  conditions  and  the  general  algorithm  for  the  solution 
cf  the  entire  boundary  value  problem.  You  are  encouraged  to  present  any  other 
Inform, atlor,  which  you  thin*  will  be  relevant. 


If  at  all  possible,  you  are  requested  to  provide  us  with  hard  copies  cf 
your  presentation  or  at  least  those  viewgraphs  responding  to  Enclosure  (1) 
approximately  one  week  before  the  meeting. 

It  is  our  intention  to  provide  those  individuals  whose  codes  are  evaluated, 
with  unique  private  Identifiers  so  that  only  the  best  codes  will  be  explicitly 
Indentlfied  during  the  comparative  evaluation.  This  Is  Intended  to  avoid 
embarrassment  to  those  Individuals  whose  codes  have  significant  deficiencies. 


We  look  forward  to  seeing  you  on  the  12th  and  13th.  If  you  have  any 
questions,  please  feel  free  to  contact  me  at  202-227-1450  vhlcn  is  my  office 
number  or  301-567-3270  at  my  heme. 


Si ncerely , 

Arthur  H.  Reed 
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WAKE-OFF  QUESTIONNAIRE 

In  order  to  complete  the  Kelvin  wake  code  evaluations,  the  committee  is  requesting 
all  participants  to  supply  answers  to  the  following  questions  regarding  the  numerical  code 
which  was  used  as  well  as  the  resulting  computations. 

COMPUTATIONAL  CODE 
Name  of  Code? 

How  many  lines  of  Coding? 

How  many  subroutines? 

What  are  the  limitations  on  the  number  of  panels  which  can  be  used? 

Describe  the  status  of  the  program  documentation. 

Does  a  user’s  manual  exist? 

What  hardware  could/has  the  Code  run  on?  With  what  modifications? 

Can  the  Wave  Resistance  Coefficient,  Cw,  be  predicted? 

Can/Does  the  Code  include  lift?  Describe  the  difficulties  of  modifying  the  Code  to 
include  lift  and/or  nonlinearities. 

Can/Does  the  Code  include  provision  for  a  propulsor/actuator  disc? 

Are  on-body  and/or  off-body  streamlines  calculated? 

What  kind  of  boundary  conditions  or  numerical  techniques  are  used  to  deal  with  the 
radiation  condition? 

Is  the  Code  suitable  for  transom-stem  type  hull  forms? 

What  experience  level  is  necessary  to  prepare  the  input  data,  run  the  program,  and 
display  the  results? 

Describe  the  ltuser-friendliness"  aspects  of  the  Code. 

KELVIN  WAVE  COMPUTATIONS 

On  what  machine  were  the  computations  run? 

How  long  (CPU  time)  did  it  take  to  run  each  case? 

How  involved,  streamlined,  and  time  consuming  is  the  input  data  preparation 
process? 

Was  a  check  done  for  convergence  during  or  after  the  computations? 

How  many  panels  were  used  (for  each  case/huli)? 

What  is  the  nature  of  the  panel  shape  and  distribution,  and  singularity'  distributions? 
If  Cu.  was  calculated,  please  provide  the  predictions  for  each  case. 
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Experimental  wave  spectrum  (or  Model  5415  at  Fn  =  0.25. 
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0.2755  SPEED  (FPS)  =  6  77  PROBE  No. 
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EXPERIMENTAL  FREE 


Fig.  B.6.  Experimental  wave  spectrum  for  Model  54 1 5  at  Fn  -  0.2755. 
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NONDIMENSIONAL  LONGITUDINAL  DISTANCE,  xKo 


1 0. 1 7  PROBE  No.  t  YFT  =  6.08 


Experimental  wave  spectrum  for  Model  5415  at  Fn  =  0.4136. 
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FARWAV  PREDICTIONS 


WAKE-OFF  QUESTIONNAIRE 
COMPUTATIONAL  CODE 

Name  of  Code?  We  have  Two  Codes,  Named  EPPACl  and  FARWAV. 

EPPAC1  provides  an  approximate  solution,  defined  explicitly  in  terms  of  the  hull 
shape  and  the  value  of  the  Froude  number,  to  the  Neumann-Kelvin  problem.  This  code 
predicts  the  flow,  both  on  the  mean  wetted  hull  surface  and  at  any  presr  ribed  points  (in¬ 
cluding  points  in  the  far  field)  outside  the  ship,  according  to  the  so-called  first-order 
slender-ship  approximation  presented  by  Noblesse  In  the  Journal  of  Slap  Research,  Vol, 
27,  pp.  13—33  (1983),  EPPACl  can  either  be  used  alone  (as  was  the  c use  for  the  calcula¬ 
tions  we  performed  for  the  DDG51  and  QUAPAW  models),  or  it  can  be  coupled  with 
(more  precisely:  its  predictions  can  be  used  as  input  to)  the  code  FARWAV. 

FARWAV  computes  the  wave  potential,  at  any  prescribed  point  behind  the  ship 
stem,  given  the  value  of  the  Froude  number,  the  hull  shape  and  the  value  of  the  velocity 
potential  on  the  mean-wetted  hull  surface.  If  the  disturbance  veloci  ty  potential  is  set 
equal  to  zero,  the  wave  potential  predicted  by  FARWAV  is  identical  to  that  predicted  by* 
EPPACl.  FARWAV  is  based  on  the  Neumann-Kelvin  theoretical  framework.  This  code 
can  be  coupled  with  ANY  code,  not  just  EPPACl,  that  predicts  the  potential  on  the 
mean-wetted  hull  surface. 

How  Many  Lint's  of  Coding  ?  How  Many  Subroutines ? 

EPPACl  and  FARWAV  each  consist  of  35  subroutines  totaling  nearly  14,000  lines,  a 
significant  portion  of  which  are  comments,  Of  these  35  subroutines  only  18,  totaling 
9.722  lines  are  actually  responsible  for  the  calculations.  The  remaining  17  subroutines  of 
4,275  lines  are  nearly  evenly  divided  between  pre-processing  and  post-processing.  The 
pre-processing  routines  manipulate  the  information  about  the  hull  into  the  format  re¬ 
quired  by  EPPACl.  The  post- processing  routines  combine  the  individual  components  of 
the  potential  and  output  than  in  a  format  ready  for  f  lotting. 

Limitations  on  the  Number  of  Panels  which  Can  be  Used? 

Dimensions  for  the  hull  panels  are  currently  set  at  150  x  21,  but  these  can  be  modi¬ 
fied  at  will  as  there  are  NO  limitations  on  the  number  of  panels  which  can  be  used  in 
these  two  codes.  The  current  dimensions  appear  to  be  sufficient  for  most  applications, 
based  on  the  convergence  tests  we  have  performed. 

Program  Documentation  ? 

The  theory  underlying  the  cedes  has  been  published,  We  are  currently  in  the  pro¬ 
cess  of  completing  some  deuile  A  numerical  studies  designed  to  refine  and  finalize  several 
numerical  aspects  of  the  codes.  Most  of  these  numerical  studies  will  be  completed  in  the 
very  near  future.  We  then  intend  to  prepare  a  series  of  DTRC  reports  presenting  a  de¬ 
tailed  account  of  the  numerical  aspects  of  the  codes  (e.g.,  methods  used  for  evaluating  the 
Green  function  and  for  integrating  it  over  a  panel)  and  comparisons  of  numerical  results 
and  experimental  data.  We  expect  that  most  of  this  work  can  be  performed  in  1988, 
although  some  particular  points  will  still  require  further  work;  and  that  will  be  presented 
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whit  completed,  We  alio  expect  tbit  this  detailed  dGeumematien  of  the  numerical  as= 
pec ts  of  the  eodii  eat  ultimately  be  published,  at  *  Conference  or  in  the  Journal  of  Ship  - 
Research, 

User  Manual? 

No  user  manual  exist!  at  this  time,  and  we  do  not  envision  preparing  one  in  the  near 
future,  The  reason  is  that  a  number  of  modifications  and  refinements  can  and  must  be 
done,  We  therefore  believe  It  would  be  most  appropriate  to  complete  the  previously 
mentioned  numerical  studies  and  documentation  reports  first,  and  then  to  modify  the 
codes  in  accordance  with  the  results  of  the  aforementioned  studies  and  with  the  expert* 
enee  of  using  the  codes,  and  finally  to  prepare  a  user's  manual, 

Hardware ? 

The  codes  have  been  run  on  APOLLO  and  MICROVAX  micro-computers,  They 
could  presumably  be  run  with  very  little  modification  on  any  virtual  memory  computer, 
Use  on  a  computer  with  limited  memory  would  require  writing  some  variables  to  scratch 
files. 

Wave  Resistance? 

The  wave  resistance  coefficient  can  readily  be  evaluated  in  both  EPPACl  and  FAR- 
WAV  since  these  codes  compute  the  wave  amplitude  function  (also  known  as  wove 
spectrum  or  Kochin  function).  A  comparison  between  the  wave  resistance  coefficients 
predicted  by  EPPACl  and  FARWAV  and  available  experimental  data  for  the  Wlgley  hull 
held  in  fixed  position  (no  sinkage  or  trim  permitted)  is  attached  to  this  questionnaire, 

The  dashed  line  corresponds  to  the  zeroth-order  approximation  to  the  wave-resistance, 
which  is  predicted  by  EPPACl.  The  solid  line  was  predicted  by  FARWAV,  v  ith  EPPACl 
used  as  input,  as  thus  corresponds  to  the  first-order  approximation.  It  is  interesting  th'.t 
the  amplitude  of  the  humps  and  hollows  in  the  resistance  curve  are  considerably  atte¬ 
nuated  in  the  Erst  approximation,  in  comparison  to  the  zeroth  approximation. 

Lift,  Nonlinearities,  Propulsor  Actuator  Disk? 

The  codes  include  neither  lift  nor  nonUneurities.  Flow  past  lifting  surfaces  can  of 
course  be  modeled  with  the  Neumann-Kelvin  theoretical  framework,  by  including  nor* 
mal-dipole  sheets  and  using  the  Kutta  condition,  This  extension  would  entail  substantial 
modifications  of  the  codes,  however, 

Nonlincarities  could  be  included  in  FARWAV,  not  in  EPPACl,  The  nonlinear  terms 
in  the  free-surface  boundary  condition  are  modeled  via  an  integral  on  the  mean-free  sur¬ 
face  plane.  The  density  of  the  free-surface  source  distribution  representing  the  nonlinear 
effects  could  be  evaluated  by  usi^g  EPPACl.  or  some  other  ruar-field  code.  This  ap¬ 
proach  could  probably  be  implemented  relatively  easily,  and  would  represent  a  fairly 
simple  extension  of  the  existing  codes,  However,  extensive  testing  would  probably  be 
required  to  ascertain  the  accuracy  of  the  calculations. 

Provision  for  a  propulsor/actuator  disk  could  probably  be  Included  in  the  codes,  at 
least  in  some  approximate  manner,  by  including  a  normal-dipole  sheet. 
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Streamlines? 

Streamlines  arc  not  calculated  in  the  current  version  of  EPPACL  These  calculations 
could  however  be  included  since  Vq  can  be  evaluated  everywhere. 

Radiation  Condition? 

The  radiation  condition  is  satisfied  automatically  and  exactly  in  these  codes  by 
means  of  the  Green  function,  which  only  permits  waves  behind  the  singularity. 

Transom-Stem  Type  Hull? 

■  The  codes  can  be  used  for  transom-stern  type  hulls,  by  assuming  the  hull  is  ex¬ 
tended  smoothly  into  a  closed  wake,  and  by  creating  the  ffaasom-stem  hull  and  its  wake 
as  a  (modified)  closed  hull.  To  be  sure,  this  approach  is  merely  approximate  and  indeed 
probably  fairly  crude.  A  more  refined  potential-flow  model  for  treating  transom-stern 
ships  should  be  formulated.  We  suspect  that  this  would  lead  to  modifications  that  might 
be  included  into  FARWAV  fairly  easily. 

Experience  Required  for  Using  the  Code?  User  Friendliness? 

The  codes  arc  pretty  straightforward  and  could  probably  be  used  by  any  engineer 
with  fairly  little  training.  Nevertheless,  the  required  input  data  could  be  simplified  by 
modifying  the  codes  somewhat.  We  would  like  to  perform  these  modifications  at  some 
later  dme. 

KELVIN  WAKE  COMPUTATIONS 
Computer? 

The  reported  calculations  for  the  DDG51  and  the  QUAPAW  models  were  per- 
formed  on  an  a\POLLO  microcomputer,  by  using  the  code  EPPAC1.  The  code  FARWAV 
was  not  used  for  these  calculations,  because  we  were  quite  pressed  to  meet  the  deadline 
and  had  not  had  any  time  for  validating  FARWAV,  Both  EPPAC1  and  FARWAV  have 
since  bien  run  on  a  MICRO  VAX  computer  for  another  hull  form  (a  Scries  60  hull). 

Computing  Time? 

The  computing  time  used  for  computing  the  velocity  potential  on  the  hull  surface 
and  at  the  free  surface  with  EPPAC1  on  an  APOLLO  workstation  was  approximately  1 
hour  per  Froude  number.  However,  we  regard  this  computing  time  as  very  excessive,  and 
we  have  now  devised  improvements  in  several  basic  aspects  of  the  code  which  we  expect 
will  reduce  the  computing  time  very  significantly  by  a  factor  of  four  or  better.  Significant 
further  improvements  can  still  be  made,  We  believe  that  we  can  have  EPPAC1  and  FAR¬ 
WAV  run  ^terac lively  on  a  microcomputer. 

Input  Data  Preparation  ? 

The  codes  require  that  the  hull  be  approximated  by  means  of  a  set  of  triangular  pan¬ 
els,  which  is  a  fairly  routine  (although  not  trivial)  task. 

Convergence? 

Checks  for  convergence  were  NOT  directly  performed  for  the  DDG51  and  the 
QUAPAW  models,  because  of  banted  time.  However,  checks  have  been  performed  for 
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the  Wigley  hull  and  a  Series  60  hull.  These  numerical  studies  on  the  effect  of  panel  size  - 
(and  other  aspects  influencing  the  accuracy  of  the  calculations)  led  to  conclusions  and 
guidelines  which  were  used  in  our  calculations  for  the  DDG51  and  QUAPAW  models. 

Panelization? 

EPPAC1  requires  only  a  relatively  small  number  of  panels;  typically  30  to  70 
“framelines"  by  15  to  18  “waterlines,"  depending  on  the  Froude  number  and  the  hull 
form.  FARWAV  requires  more  panels,  because  the  velocity  potential  varies  more  rapidly 
then  the  hull  form  (especially  at  low  Froude  numbers).  However,  we  recently  obtained  a 
modified  mathematical  expression  for  evaluating  the  wave-amplitude  function  (Kochin 
function)  which  requires  only  the  tangential  velocity  components  (not  the  velocity  poten¬ 
tial)  on  the  mean  wetted  hull  surface  as  input  (as  well  as  the  hull  shape  and  the  Froude 
number,  of  course).  This  fact,  together  with  mathematical  properties  of  the  modified  ex¬ 
pression,  lead  us  to  expect  that  fewer  panels  will  be  required. 

The  hull  is  approximated  by  means  of  flat  triangular  panels,  within  which  the  veloc¬ 
ity  potential  is  assumed  to  vary  linearly  (in  the  modified  version  of  FARWAV  based  upon 
the  modified  expression  for  the  Kochin  function  noted  previously,  linear  variation  will  be 
assumed  lor  the  tangential  velocity  components). 

Wave  Resistance? 

The  wave  resistance  of  the  DDG51  and  QUAPAW  models  was  not  computed,  be¬ 
cause  it  was  not  required.  This  information  could  readily  have  been  obtained  since  the 
wave-spectrum  function  was  determined. 


MODEL  d5415 


FARWAV  prediction  of  wave  contours  for  Model  5415  at  Fn  =  0.25. 


MODEL  d5415 
Fn  =  0.2755 
A  {T]  Ko)  =  O.Ot 


FARWAV  prediction  of  wave  contours  for  Model  541 5  at  Fn  =  0.2755. 


FARWAV  prediction  of  wave  cut  for  Model  541 5  at  Fn  =  0.2755. 
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Fig.  C.6.  FAR  WAV  prediction  of  wave  spectrum  for  Model  5415  at  Fn  =  0.2755. 


MODEL  d54 15 
Fn  =  0.4136 


F  ARWAV  prediction  of  wave  contours  for  Model  541 5  at  Fn  =  0.4136. 


0.41 3S 


Fig.  C.8.  FARWAV  prediction  of  wave  cut  for  Model  541 5  at  Fn  =  0.4136. 


MODEL  3531 
Fn  =  0.2131 
A  (T\  K0)  =  0.03 


FARWAV  prediction  of  wave  contours  for  QUAPAW  at  F(I  =  0.2131 


G.2T3I 


FAR  WAV  precfictian  of  wave  cut  for  QUAPAW 


MODEL  353t  F„  =  0.213t 


FARWAV  prediction  of  wave  spectrum  for  QUAPAW  af  Fn  -  0.2 T 31! 


MODEL  3531 


FARWAV  prediction  of  wave  contours  for  QUAPAW 


0.2500 
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Fig.  C.14.  FARWAV  prediction  of  wave  cut  for  QUAPAW  at  Fn  =  0.250 


MODEL 


FARWAV  prediction  of  wave  spectrum  for  QUAPAW  al 


MODEL  3531 


F  ARWAV  prediction  of  wave  contours  for  QUAPAW  at  Fn  -  0.3197. 


2.0- 


MODEL  3531  Fn  =  0.3197 


FAFtWAV  prediction  of  wave  spectrum  for  Q  JAPAW  at  Fn  =  0.3 1 97. 


*ZiZ =1-. — r 
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APPENDIX  D 
NK-2  PREDICTIONS 


PANELIZED 

GEOMETRY 


I 
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COMPUTATIONAL  COPES  _ 


CODE  NAME: 

NK_:2 

NEAR -WAKE 

FAR-WAKE 

Lines: 

2350 

(+  IMSL  EQ.  Solver) 

1600 

800 

Subroutines: 

46 

23 

12 

Number  of  Panels: 

Unlimited 

Unlimited 

Unlimited 

Documentation: 

Poor 

Poor 

Poor 

User's  Manual: 

No 

No 

No 

’  udware: 

VAX  760  &  8650 
CRAY  X-MP 

VAX  750  &  8650 
CRAY  X~MP 

VAX  750  &  8650 
CRAY  X-MP 

Cw: 

No 

No 

Yes 

Lift: 

Yes 

Yes 

Yes 

Actuator  Disc: 

No 

No 

No 

Streamlines: 

Nr 

Yes 

No 

Rauiauon 

Condition: 

Havelock  Source 

Haveiock  Source 

Havelock  Source 

Transom  Steam: 

n 

Expenence 

Level: 

Moderate 

High 

Moderate 

User  Friendliness: 

Moderate 

Low 

Moderate 
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KELYTN  WAVE  COMPUTATIONS 


Code  Name:  NK-2  NEAR-WAKE 

Machine:  VAX  86SC  VAX  8650 

CPU  Cost:  $52/hour  (Night  Batch  Rate) 

Seconds  of  CPU  Time  (Cost): 

QUAPAW  (432  Panels) 


PAR-WAKE 
VAX  8650 


10.0  kts 

11.7  kts 

15.0  kts 

1150  (SI  7) 

1200  ($17) 

1300  ($19) 

4200  (S61) 

1700  (S25) 

1600  ($23) 

ICO  (SI) 
100  (SI) 
100  ($1) 

Model  5415(376  Panels) 

18.1  kts 

20.0  kts 

30.0  kts 

1300  ($19) 

1400  ($20) 

2700  ($39) 

4400 ($64) 

4400  ($64) 

4400  (S64) 

100  ($1) 
100 ($1) 
100  ($1) 

Input  data  Preparation; 

GEOPAN  (SAIC  panelization  code) 

Convergence:  Will  be  discussed 


Panels:  Will  be  discussed 


En 

Cw 

QUAPAW 

.213 

1.37 

S-9000 

.250 

6.40 

.320 

7.92 

Model  5415 

.250 

0.51 

5^32200 

.276 

1.39 

.414 

4.39 

68 


CONTOUR  INTERVAL  =  0.020  x  1/K0 


WAVE  ELEVATION 


WAVE  ELEVATION 


Xx  WAVE  NUMBER 


WAVE  ELEVATION 


WAVE  ELEVATION 


WAVE  ELEVAT ION 


-2  prediction  of  wave  cut  for  Model  541 5  at  Fn  =  0.2755 


FREE-WAVE  SPECTRUM  (SHARMA) 
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WAVE  ELEVATION 


77 


prediction  of  wave  contours  [-0.5  <  2wLWL  <  4)  for  Model  5415  at  F„  =  0.4136. 


CONTOUR  INTERVAL  =  0.020-  t/K 


73 


E  SPECTRUM  (SHARMA) 


spectrum  for  Model  5415  at  Fn  =  0.4136. 


WAVE  ELEVATiON 


-A1 


SPECTRUM  (SHARMA) 


WAVE  ELEVATION 


WAVE  ELEVATION 


WAVE  ELEVATION 


prediction  of  wave  cut  for  QUAPAW 


E-WAVE  SPECTRUM  (SHARMA) 


aonindwv  ivnoisn3^iq-non 
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0  0 


WAVE  ELEVATION 


CONTOUR  INTERVAL  =  0.020 


WAVE  ELEVATION 


prediction  of  wave  cut  for  GUAPAV 


SPECTRUM  (SHARMA) 

Pxk2  +  Gxh2 


spectrum  fur  QUAPAW  at  Fn  =  0.3197. 


APPENDIX  E 

DOCTORS  N-K  PREDICTIONS 


Program  -  DOCTORS.N-K 
A  Neumann-Kelvln  code  using 

the  source  distribution  technique  KELVIN  WAVE  COMPUTATIONS 


9' 


no  limitation  on  number  of  field  points 


WAVE  CUT 


Fig.  E.3.  DOCTORS  N-K  prediction  of  wave  contours  for  QUAPAW 


MODEL  353 f  Fn  =  0.2131 


Fig.  E.  4.  DOCTORS  N-K  prediction  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.21311 


2.0 


MODEL  3531 


MODEL  353t 


FJg.  E.7.  DOCTORS  N-K  prediction  of  wave  contours  for  QUAPAW 


MODEL  3531  Fn  =  0.25 
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Fig.  E.fe.  DOCTORS  N-K  predictiorr  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.25. 


I 

o 


in 

o 
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WAVE  CUT 


DOCTORS  N-K  prediction  of  wave  contours  for  QUAPAW 


MODEL 


DOCTORS  N-K  prediction  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.3197. 


APPENDIX  F 


SWIFT  PREDICTIONS 
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MODEL  5415  Fa  =  0.25 


MODEL  5415 


LONGITUDINAL  DISTANCE  FROM  THE  BOW  IN  X*g/U**2 


MODEL  54T5  Fn  =  0.2755 


MODEL  5415  Fn  =  0.2755 
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MODEL  5415 


o 


3**n/6*z  Nl  d31VM  W1V0  01  3All\n3ti  1HDI3H  3AVM 
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MODEL  5415  Fn  =  0.4136 


MODEL  5415  Fn  =  0.4136 
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MODEL  5415 


S**n/6*Z  Nl  H31VM  W1V0  Ol  3AI1\T!3H  1H0I3H  3AVM 
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MODEL  3531 


MODEL  3531 


SWIFT  prediction  of  wave  cut  forQUAPAW  at  F0  =  0.25. 


0.3197 
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APPENDIX  G 
NKSHXP  PREDICTIONS 
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WAKE-OFF  QUESTIONNAIRE 
Response  from  MIT 


COMPUTATIONAL  CODE 


(1)  Name  of  Code? 

NKSHIP 

(2)  Who  many  lines  of  Coding? 

1550  plus  standard  routines  (RPAN,  GEOM,  etc.). 

(3)  How  many  subroutines? 

12  plus  standard  routines. 

(4)  What  are  the  lirr,itations  on  the  number  of  panels  which  can  be  used? 

Specified  by  available  machine  memory.  If  direct  solver  is  used  on  VAX-750  1500/3000 

(5)  Describe  the  status  of  the  probram  documentation . 

Under  development,  not  commercialized. 

(6)  Does  a  user’s  manual  exist? 

No. 

(7)  What  hardware  could/has  the  code  run  on?  With  what  modifications? 

On  any  machine  that  supports  FORTRAN  77,  without  any  modifications. 

(8)  Can  the  Wave  Resistance  Coefficient,  Cm  ,  be  predicted? 

Yes,  but  the  prediction  is  not  reliable. 

(9)  Can/Does  the  Code  include  lift?  Describe  the  difficulties  of  modifying  the  code  to 
include  lift  and/or  non-linearities. 

No,  at  least  at  this  stage. 

(10)  Can/Does  the  Code  include  provision  for  a  propulsor/ Actuator  disc? 

No,  at  least  at  this  stage. 

(11)  Are  on-body  and/or  off-body  streamlines  calculated? 

Yes,  on  the  body  and  on  the  free  surface. 

(12)  What  kind  of  boundary  conditions  or  numerical  techniques  are  used  to  deal  with  the 
radiation  condition? 

The  freesurfaee  elevation  and  slope  vanish  at  the  upstream  boundary  of  the  free  surface 
computational  domain,  while  no  conditions  are  imposed  at  the  downstream  boundary'. 
[Note:  the  shape  function  used  for  the  approximation  of  the  potential  is  a  quatratic 
spline  that  needs,  in  general,  two  end-conditions  in  each  direction.) 

(13)  Is  the  Code  suitable  for  transon-stern  type  hull  forms? 

It  has  not  been  tried  yet. 

(14)  What  experience  level  is  necessary  to  prepare  the  input  date,  run  the  program  and 
display  the  results? 

The  code  is  not  commercialized  yet. 

(15)  Describe  the  “ user-friendliness “  aspects  of  the  code. 

The  code  is  not  commercialized  yet. 
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KELVIN  WAVE  COMPUTATIONS 


(1)  On  what  machine  were  the  compulations  run? 

VAX-750  and  CRAY-XMP. 

(2)  How  long  ( CPU  time)  did  it  take  to  run  each  case? 

3  hrs  on  a  VAX-750. 

(3)  How  involved,  streamlined,  and  time  consuming  is  the  input  data  preparation  process? 
The  code  is  not  commercialized  yet. 

(4)  Was  a  check  done  for  convergence  during  or  after  the  computations? 

Yes  -  we  think  that  this  is  essential. 

(5)  How  many  panels  were  used  (for  each  ease/hull)? 

905,  of  which  585  are  on  the  free  surface  and  320  on  the  body. 

(5)  What  is  the  nature  of  the  panel  shape  and  distribution  and  singularity  distributions? 
Plane  quadrilaterals,  quadratic  spline. 

(6)  If  Cw  was  calculated,  please  provide  the  predictions  for  each  case. 

Yes,  but  without  being  reliable. 
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0.3197 


NKSHtP  prediction  of  wave  contours  forQUAPAW  at  Fn  =  C.3197. 


ON  FS 
ON  BODY 


NiKSHtP  orediction  of  wavs  profiles  along  QUAPAW 


ON!  BODY 


NKSH1P  prediction  of  wave  profiles  along  QUAPAW 


ON  BODY 


g.  G.5.  NKSHIP  prediction  of  wave  profiles  along  QUAPAW  hu!'  at  Fn  =  0.3197. 


APPENDIX  H 


FLOPAN  PREDICTIONS 
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FLOPAN  half-plane  panel  model  of  QUAPAW. 
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FLOPAN  prediction  of  wave  contours  forGUAPAW 


CONTOUR  INTERVAL 


-to 
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XxKc 


2  0 


°X  *  VIS 
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of  wave  cut  for  GUAPAW  at  Fn  =  0.31 97. 


0.4- 
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FLOPAN  prediction  of  wave  profile  along  QUAPAW 


0.4 
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FLOFAN  prediction  of  wave  profile  along  QUAPAW  hull  at  F„,  =  0.25. 


0.4- 


FJg.  H.10.  FLOPAN  prediction  of  wave  profile  along  QUAPAW  hull  at  Fn  =  0.31 97. 
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APPENDIX  I 
XYZFS  PREDICTIONS 
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XYZFS  PROGRAM 

•  GENERAL  DESCRIPTION 

•  Originally  developed  by  Charles  Dawson 

•  Further  developed  by  Bill  Cheng  and  Janet  Dean 

•  Ranking  source  panel  method 

•  Computes  3D  flow  fields  around  a  ship  in  forward 
motion 

•  2000  lines 

•  Seven  segments  plus  numerous  preprocessors  and 
postprocessors 

.  XYZFS  runs  on  CRAY  IS  and  CRAY  XMP-24. 

Number  of  panels  limited  only  by  computer  memory 
(currently  1500). 

•  Preprocessors  and  postprocessors  run  on  /polio 
workstations. 

•  COMPUTATIONAL  METHODS 

•  Double  model 

•  Free  surface  solution 

.  •  Direct  solution  method  (Gaussian  Elimination) 

•  Radiation  condition  satisfied  by  a  four-point  upstream 
finite  differrence  operator. 

•  Velocity  induced  by  a  panel  is  computed  using  exact 
integration  and  multiple  expansion  methods  developed 
by  Hess  and  Smith  (1962). 

•  USER  REQUIREMENT  /  BACKGROUND 

•  Input  preparation  —  ship  geometry  and  lines  drawing 

•  Program  execution  —  previous  computer  usage 

•  Results  display  —  computer  graphics  (2D,  3D  and 
contour  plots) 

•  Results  analysis  —  physics  of  ship  waves 

•  USER-FRIENDLINESS 

•  User’s  manual  available 

•  Provisions  for  error  checkings  of  pane)  data  such  as 
aspect  ratio,  outward  normal,  deviations  ‘roni  a  flat 
panel. 

•  RESEARCH  AREAS  OF  PANEL  METHODS 

«  Nonlinearities  —  Nonlinear  free  surface  w  r  c  has  been 
reported  by  Ogiwara  and  Maruo  (19S5),  Xia  (19S6), 
Buigr.relli  (19S7),  and  Ni  (19S7). 

•  Lilting  Surface  —  Work  on  linear  free  surface  with  lift 
has  been  reported  by  Slooff  (1984)  and  Larsson  and 
Xia  (19S6). 
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KELVIN  WAKE  COMPUTATIONS  USING.  XYZFS 

•  Computations  for  QUAPAW  and  DDGS1 

•  Machine  used:  CRAY  XMP-24 

•  CPU  time:  approximately  1  minute  per  Froude  number 

•  Number  of  panels  used: 


HULL 

FREE  SURFACE 

QUAPAW 

192 

1220 

MODEL  5415 

320 

1128 

*  Quadrilateral  panels  of  constant  source  strength  were 
used.  The  panels  are  organized  by  station  lines;  points 
on  a  station  line  are  selected  by  specified  percentages 
in  girth  using  B-splines. 

•  Convergence  check  was  made  as  the  width  and  length 
of  the  FS  computational  domain  were  changed. 
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XYZFS  MODEL 
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XYZFS  MODEL  5415  Fn  =  0.25 


Fig.  1.2.  XYZFS  prediction  of  wave  cut  for  Model  541 5  al 


MODEL 


XYZFS  prediction  of  wave  spectrum 


XYZFS  MODEL 


XYZFS  MODEL  5415  Fn=  0.2755 
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XY2FS  MODEL  5415  Fn  =  0.2755 


XYZFS  prediction  of  wave  cut  for  Model  541 5  at  Fn  =  0.2755. 


XYZFS  MODEL  5415  Fn  =  0.2755 


XYZFS  prediction  of  wave  spectrum  for  Mode!  54 1 5  at  Fn  =  0.2755. 


XYZFS  MODEL  541 5  Fn  =  0.4136 


XYZFS  MODEL  541 5  F„=  0.4136 


XYZFS  MODEL  541 5  Fn=  0.4136 


XYZFS  prediction  of  wave  cut  for  Model  541 5  at  F0  =  0.4136. 


MODEL 


XYZFS  prediction  of  wave  spectrum  for  Model  54t5  at  Fn  =  0.4136. 


XYZFS  MODEL  3531  F„=  0.2131  Ar[Ko=0.01 


162 


XYZFS  MODEL 


XYZFS  prediction  of  wave  cut  for  QUAPAW  al 


XYZFS  MODEL  3531  Fn  =  0.2131 
(WITH  TRUNCATION  CORRECTION) 


XYZFS  prediction  of  wave  spectrum  for  GUAPAW 


XYZFS  MODEL  3531 


XYZFS  MODEL 


XYZFS  MODEL  353 


XYZFS  prediction  of  wave  cut  for  QUA.PAW 


XYZFS  MODEL 
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Fig.  1.15.  XYZFS  prediction  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.25. 


Q-3t97 


XYZFS  MODEL  3531  Fn=  0.3197 


0.3197 
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APPENDIX  J 

PROPRIETARY  CODE  PREDICTIONS 


I 
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PROPRIETARY  COMPUTER  CODE 


*  SOLUTION  PROCEDURE 

*  DOUBLE  MODEL  SOLUTION 

.  CREATE  FREE  SURFACE  GRID 

•  FREE  SURFACE  SOLUTION 

ITERATIVE  MATRIX  SOLVER 


•  FLAT  QUADRILATERAL  PANELS 

•  HULL  -  CONSTANT  RANKINE  SOURCES 

•  FREE  SURFACE 

CONSTANT  RANKINE  SOURCES 
CONSTANT  DOUBLETS 

•  RADIATION  CONDITION 

UPSTREAM  FINITE  DIFFERENCE  OPERATOR 

•  CAPABILITIES 

.  COMPUTES  WAVE  RESISTANCE 

*  INCLUDES  LIFT 

.  SOME  ERROR  CHECKING 
.  NO  SINKAGE  /  TRIM 
.  NO  STREAMLINE  OPTION 
.  NO  PROPULSOR  /  ACTUATOR  DISC 
.  MAY  BE  SUITABLE  FOR  TRANSOM  STERNS 
NOT  A  DOCUMENTED  OPTION 

•  CODING  INFORMATION 

•  FORTRAN 

.  15,000  LINES 
.  80  SUBROUTINES 
.  VAX  OR  CDC  COMPUTER 

•  DOCUMENTATION 

.  SEVERAL  REPORTS  ON  THE  METHOD 

♦  USER’S  MANUAL 

.  CODE  SPARSELY  COMMENTED 
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Proprietary  Code  prediction  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.2131 


182 


'e  cut  for  QUAPAW  at  Fn  =  0.25. 


Proprietary  Code  prediction  of  wave  spectrum  for  QUAPAW 


QUAPAW  a!  Fn  =  0.3197. 


Proprietary  Code  prediction  of  wave  spectrum  tor  QUAPAW  at  Fn  =  0.31 37. 


APPENDIX  K 
SWIMFS  PREDICTION 
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SWIM/SWIMFS  COMPUTER  CODE  —  SWIM  SWIM/SWIMFS  —  DOCUMENTATION 


5  *  S 

■§  S  S’ 


m 

*1  5 

*5 

J  w  a 


553  lines  of  FORTRAN  (excluding  commeuls) 


SWIM/SW1MFS  —  USER-FRIENDLINESS 


•  No  checking  of  geometric  input  — »  prints  centroid,  normal, 
and  area  of  each  panel. 

«  No  informative  messages  if  something  goes  wrong  or  is  likely 
to  go  wrong  in  the  computation. 

•  To  prepare  input  — 

•  To  run  the  program  — 

•  To  display  results  — 


SWIM/SWIMFS  —  computations 

*  Sunk/trimmed  Quapaw  data  from  XYZFS 

•  CPU  time:  about  50  minutes  per  spee,d  (CR\Y-XMP) 

*  152  (S  X  24)  panels  on  half  the  hull 

•  1200  to  1250  free-surface  points  for  contour  plotting  (only 
inside  the  19"  cusp  lines)  —  on  the  Quapaw: 


Speed  (knots) 

pis. /wavelength 

10.0 

13 

11.7 

S 

15.0 

6 

•  Rectangular  panels.  Singularity  distribution  cf  Havelock 
sources  on  panel  collapsed  onto  line  segment  through 
centroid  of  panel. 

•  No  convergence  check. 

•  Wave  spectra  —  from  (denser)  wave  cut  data 
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MODEL  353' 
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MODEL  3531  F„  =  0.2131 


SWIMFS  prediction  of  wave  cut  far  QUAPAW 


MODEL  3531  Fn  =  0.25  AfTtKJ  -  0.02 


MODEL  353 T  Fn  =  0.25  AfrtKJ  =  0.02 


MODEL 


SWIMFS  prediction  of  wave  cut  forQUAPAW 


MODEL  3531  Fn  =  0.3197  A|r{K0]  =  0.05 


MODEL  353t  Fn  =  0.3197  AfaKo]  =  0.05 


MODEL  353 1  Fn  =  0.3197 


SWtMFS  prediction  of  wave  cut  for  QUAPAW  at  Fn  =  0.3197. 


MODEL  3531  Fn  =  0.3197 
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Fig.  K.9.  SWIMFS  prediction  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.3197. 
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APPENDIX  L 

SAIC  SLENDER-SfflP  SUBMITTED  AFTER  WAKE-OFT 


WAVE  ELEVATION 
CONTOUR  INTERVAL  =  0. 


sfender-ship  prediction  of  wave  contours  (-0.5  <  2x/LWL<  4} 


WAVE  ELEVATION 
CONTOUR  INTERVAL  -  0.020  x 


IC  slender- ship  prediction  of  wave  contours  (4  <  2x/LWL<8) 


WAVE  ELEVATION 
VELOCITY  =  T8.1  kny/lV 


(MODEL  54t5) 

3.1  kn 
>25 

P  APPROXIMAT  (ON 


FREE-WAVE  SPECTRUM  (SHARMA) 
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Fig.  L.4.  Post  wake-off  SAIC  slender-ship  prediction  of  wave  spectrum 


WAVE  ELEVATION 
CONTOUR  INTERVAL  =  0.1 


lode!  5415  at  Fn  =  0.2755. 


WAVE  ELEVATION 
CONTOUR  INTERVAL  =  0  020  x  l/K 
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for  Mode!  54  f  5  at  Fn  =  0.2755. 


WAVE  ELEVATION 
VELOCITY  =  18.1  kn  y/L  WL 
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WAVE  ELEVATION 
CONTOUR  INTERVAL  =  0.020  x  t/K 


(S' 


215 


ELOCITY 


E  SPECTRUM  (SHARMA) 


WAVE  EL  EVATION 
CONTOUR  INTERVAL  =  0.020  x  t/K 


QUAPAW  at  Fn  =  0.2131. 


WAVE  ELEVATION 
CONTOUR  INTERVAL  =  0  020  x  t/K 
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Fig.  L.14.  Post  wake-off  SAIC  siender-ship  prediction  of  wave  contours  (4  <  2x/LWL  <  8) 
for  QU APAW  at  Fn  =  0.21 3  f . 
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Fig.  L.17.  Post  wake-off  SAIC  s!ender-ship  prediction  of  wave  contours  (-0.5  <  2x/LWL  <  4) 
for  QUAPAW  at  Fn  =  0.25. 


WAVE  ELEVATION 
CONTOUR  INTERVAL  =  0.020  x  1/K 


Post  wake-off  SAIC  slender-shfp  prediction  of  wave  contours  (4  <  2x/LWL  <  8) 
for  GUAPAW  at  Fn  =  0.25. 


WAVE  ELEVATION 
VELOCfTY  =  t1.7kn  y/LWL  = 


sfender-ship  prediction  of  wave  cut  for  QUAPAW 


E  SPECTRUM:  (SHARMA) 


WAVE  ELEVATION 


st  wake-off  SAIC  slender-ship  prediction  of  wave  contours  (-0.5  <  2x/LWL  <  4) 
QUAPAW  at  Fn  =  0.3197. 


WAVE  ELEVATION 
CONTOUR  tNTERVAL  =  0 


QUAPAW  at  Fn  =  0.31 97. 


WAVE  ELEVATION 
OCITY  =  15.0  kn  y/LWL 


-REE-WAVE  SPECTRUM  (SHARMA) 


Summary  of  Kelvin  Wake  computations  performed  by  SAIL  in 
conjunction  with  the  wake-off. 
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APPENDIX  M 


FLOPAN  PREDICTIONS  SUBMITTED  AFTER  WAKE-OFF 
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2.0- 


234 


236 


APPENDIX  N 

XYZFS  PREDICTIONS  SUBMITTED 
AFTER  WAKE-OFF 
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XYZFS  prediction  of  wave  spectrum  for  QUAPAW  at  Fn  =  0.2131 


prediction  of  wave  contours  for  QUAPAW 


XYZFS  prediction  of  wave  spectrum  for  QUAPAW 
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